Photonic spectrum of bichromatic optical lattices 
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We study the photonic spectrum of a one-dimensional optical lattice possessing a double primitive 
cell, when the atoms are well localized at the lattice minima. While a one-dimensional lattice with 
a simple Wigner-Seitz cell always possesses a photonic bandgap at the atomic resonance, in this 
configuration the photonic transmission spectrum may exhibit none, double or multiple photonic 
bandgaps depending on the ratio between the interparticle distance g inside the cell and the cell size 
a. The transmission spectra of a weak incident probe are evaluated when the atoms are trapped in 
free space and inside an optical resonator for realistic experimental parameters. 



I. INTRODUCTION 

Ultracold atoms in optical lattices constitute a 
paradigmatic system, which allows one control over sev- 
eral parameters, thereby mimicking dynamics typical of 
condensed matter systems [HQ- A remarkable feature 
of optical lattices is that the bulk periodicity is here con- 
trolled by engineering the geometry of the propagating 
beams, which determine the light potentials. Differing 
from ordinary crystals in condensed matter, the size of 
the Wigner-Seitz cell is of the order of the light wave- 
length. One consequence is that the light, coupling with 
the atomic transitions, is also diffracted by the crystal- 
lyne structure which the atoms form [3[. 

It has been observed that the modulation of the atomic 
density in these systems, and hence of the refractive in- 
dex, makes optical lattices a photonic bandgap mate- 
rial [i[ . Theoretical works studied the photonic bandgap 
for one-dimensional and three-dimensional atomic struc- 
tures [1, [H, H, 0, H, El • Besides potential applications for 
nonlinear optics, the full understanding of these proper- 
ties is important in order to measure the quantum state 
of cold atoms using light thereby opening interest- 
ing avenues in quantum information processing for imple- 
menting photonic interfaces using ultracold atoms [111 ]. 

In this work we study theoretically the photonic prop- 
erties of biperiodic optical lattices, in a setup similar to 
the ones realized experimentally in [l2|, 0, [H| ■ We fo- 
cus on a one-dimensional configuration, and develop a 
full quantum model for the light-matter interactions, as- 
suming that the atoms are well localized at the lattice 
minima. The photonic spectra and the probe transmis- 
sion are evaluated when the optical lattice is in free space 
and inside a standing wave optical resonator, as a func- 
tion of the interparticle distance g inside the primitve 
cell. 

This article is organized as follows. In Sec. |TT]the the- 
oretical model is described and the basic approximations 
are introduced. In Sec. IIIII the photonic spectra are re- 
ported and discussed, and the transmission spectra for 
a weak probe are evaluated in Sec. HVl In Sec. [V] the 
photonic properties of bichromatic optical lattices inside 
a resonator are analyzed. The conclusions are presented 



in Sec. ED 



II. THEORETICAL MODEL 

The physical system we consider is a one-dimensional 
periodic distribution of atoms in a light potential with a 
double primitive cell. We assume a sequence of N atoms 
of mass m in a standing wave created with lasers along 
the x direction. We denote the atomic positions along x 
by Xj, with j = 1, . . . ,N. Denoting by a the size of the 
Wigner-Seitz cell, the positions are given by 

Xj = la for j = 2£ , 

= £a + g for j = 21+1, 

where I = 0, 1, . . . , M — 1, and M = N/2 is the number 
of cells (assuming N even for convenience). In the case 
here discussed, we set the size a = A, where A is the 
wavelength of the light which interacts with the dipolar 
transitions of the atoms. 

Such configuration can be experimentally realized by 
using a monochromatic standing wave with wavelength 
A, to which two laser beams are superposed, such that 
they are rotated by angles of 60° and 120° with respect to 
the axis of the lattice, as shown in Fig. Q] Upon setting 
the relative phases, the resulting potential for the atoms 
has the form 



U(x) oc 1 cos 2 (fcr/2) + cos 2 (fca;) 



(1) 



and the distances between adjacent wells are c?i = A(l — 

1 — fl 2 \ \ — B 2 

— acos —§-) and d>2 = — acos —j-, with d\ + d 2 = A. An- 
other possible realization is found by superposing two 
laser beams along the x-axis, with a half frequency [l5j 
or with a three fourth frequency respect to the fre- 
quency of the main lattice [l6| . Upon setting the relative 
phases, the four-atomic elementary cell has the structure 
rfi — c?i — g?2 — ^2 (with di + d.2 = A) and the crystal has es- 
sentially the same spectral properties than the biatomic 
one considered in this paper. In this work we will also 
consider biperiodic lattices, where the atoms composing 
the Wigner-Seitz cell may have different scattering prop- 
erties, for instance, they can belong to different species 
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FIG. 1: A possible optical realization of the double period 
ID lattice, considered in this work, can be obtained by using 
a monochromatic standing wave with wavelength A, to which 
two laser beams are superposed, such that they are rotated by 
angles of 60° and 120° with respect to the axis of the lattice. 
Here, di = g, d± + — X and the size of the Wigner-Seitz 
cell is a — X. 



or belong to the same species but are prepared in dif- 
ferent hyperfine states. Under the assumption that the 
frequencies of the two transitions are sufficiently close 
to allow significant coupling with the same probe, such 
lattices could be realized with linearly polarized coun- 
terpropagating beams, controlling the angle between the 
polarization [171 ] . 

In developing the theoretical model we will make the 
following assumptions : (i) The atoms are well localized at 
the lattice points, and the size of the atomic wave packet 
is very small with respect to the laser wavelength (Lamb- 
Dicke regime [Hj]). This situation can be realized when 
the atoms are deep in the Mott-insulator quantum state 
of the biperiodic potential [l| . We will treat the atoms as 
pointlike, hence considering the response of the medium 
at lowest order in the expansion of the size of the wave 
packet over the wavelength, (ii) We consider the coupling 
of the lattice only with modes of the electromagnetic field 
which propagate along the lattice. This approximation 
is valid when the atoms are placed, for instance, inside a 
bad cavity with sufficiently large cooperativity [l9| or a 
hollow-core fiber [13,[2l|. 

A. Hamiltonian 

The Hamiltonian H for the total system, composed by 
the atomic spins aligned along the x axis, and the modes 
of the electromagnetic field propagating along x, reads 

H = -ffdip + -ffcmf + #int ■ 

Here, 

-ffdip = E hujO^Uj (2) 

3 

describes the array of dipoles, with j = 1, . . . , N label- 
ing the atoms, where a 3 is the dipolar lowering operator 
and crt its adjoint, whereby the relevant transitions of 



the atoms at the even (odd) sites, j = 21 (j = 21 + 1), 
have dipole moments T>\ (D2) and transition frequency 
u>i (u>2). The Hamiltonian for the modes of the electro- 
magnetic field takes the form 

= J2 E ^4 n)t 4 n) , (3) 

k n=l,2 

where the operators and annihilate and create 
a photon in the electromagnetic field mode at frequency 
u>k, wave- vector k — ujk/c, with c light velocity, and po- 
larization e„ = £1,62, with ei 1 £2 1 x. Finally, the 
interaction between photons and atoms is described by 
the Hamiltonian term (in Coulomb gauge) 

H int = - e y2^-A ± ( Xj ), (4) 
3 J 

where m 3 is the mass of the atom j, pj is the momentum 
operator of the electron at atom j, 

TTtj / t \ 

/'/ i—~jDjUt, g s \. (5) 
and A\_ (r) is the transverse vector potential 

^M=EE&^/^4 (6) 

k n— 1,2 

with e the vacuum permittivity and V the quantization 
volume. Note that we have used a plane waves decom- 
position, applying periodic boundary conditions at the 
lattice borders. 



B. Weak excitation regime 

In this work we consider that the atomic transitions 
are driven well below saturation, and correspondingly the 
mean number of photonic excitations inside the system 
is much smaller than the total number of spins N. In 
this regime we use the Holstein-Primakoff representation 
of spin operators [2^], and expand all operators at the 
lowest orders in the powers of bosonic operators bj , 

a] = b){l-b)b^^b)(l-b)b 3 /2), (7) 

a- = (1 - b]b 3 )^ b 3 ~ (l - b 3 , (8) 

*3=-l+bfr- (9) 

In this representation, the Hamiltonian for the dipoles 
becomes the sum of N harmonic oscillators, 

j 

where we discarded the constant term. The interaction 
term reads (in the Rotating Wave Approximation) 

H int =HU+H®, (11) 
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with 

= E n ^k b ) a k ] e ikXj + H - c - . (12) 

j,k,n 

while describes the corrections beyond the linear 

response. Here, Q^jj is the coupling strength of the atom 
j with the mode (k, n) and is given by 

We will consider the limit in which we can truncate the 
expansion and approximate H lnt ~ thereby re- 

stricting to the case in which the medium polarization 
is linear in the electric field amplitude. 



where the quasi-momentum verifies the relation k = q + 
G. In this form, the Hamiltonian can be rewritten as the 
sum of M Hamiltonian terms, H — J2 q eBZ Hqi wnere 



H q = fiuib\b q + Tiuj 2 d\d q + Uq+Ga q +G a q U +G 

G,n 

G,n 

(18) 

This separation is only valid in the linear regime since 
saturation effects, described by mix the manifolds 

identified by the Hamiltonian terms H q . 



C. Spin waves 



Given the periodic structure, it is convenient to de- 
scribe the dipolar excitations in momentum space. At 
this purpose, for a sufficiently large crystal we assume 
Born-von Karman periodic boundary conditions, and 
consider the spin- wave excitations 

M-l 

M-l 

d q = 7=6-^ £ b 2i+1 c-^« , (15) 

with q the wave vector sweeping the first Brillouin zone 
(BZ). We denote by 

Go = 2ir/a 

the elementary vector of the reciprocal lattice, such that 
the interval of the first BZ is [—Go/2, Go/2]. Using the 

relation X^o* ex P(K9 — Q r )^ a ) = M5 qq i where the equal- 
ity q = q' is defined modulus a vector G of the reciprocal 
lattice, the Hamiltonian terms transform as 

#di P = h{uib\b q +uj 2 d\d q ) , (16) 

qeBZ 

^ (1) -E E ^{^ +G bl + e^gS +G d q )a q % 

G,n qeBZ 

+H.c, 

(17) 



III. THE PHOTONIC BAND STRUCTURE 



In this section we study the photonic spectrum of the 
biperiodic structure assuming that the polarization of the 
incident light, say ei, is parallel to the dipole moments 
T>i and T> 2 . Hence, we drop the polarization superscripts 
where they appear. The photonic band structure is found 
by solving the Hcisenberg equations of motion for each 
Hamiltonian block H q , 



dq+G — —lU q+ G(lq+G 

-iVM {Gl q+G b q - o-' lGe Gt q+G d q ) , (19) 
6g = -iwi6 g -i\/M^0i i9+G ag + G, (20) 

G 

d q = -iu 2 d q - iVM J2 j Ge G2,q+Ga q+G , (21) 

G 



where u q = c\q\. Hence, the spin wave at wavevector q 
couples in principle with all photonic modes at wavevec- 
tors q + G. Nevertheless, only the coupling of the atomic 
transition with the quasi-resonant modes at wavevectors 
Q = ±Go is significant. Taking into account only the rel- 
evant coupling, we can solve analytically the eigenvalue 
problem around q ~ in the BZ and assuming u)\ = uj 2 . 
In this limit one finds four eigenfrequencies, 



± 



i 



LUQ — LOl 



+ MQ 2 



2|{?1,Q<?2,<2| 



sin 2 GqQ 



(22) 
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where Q = ^/\Gi,q\ 2 + \G2,q\ 2 and j = 1, 2. They deter- 
mine the edges of two photonic bandgaps, one at the fre- 
quencies between v\ t - and V2,- and the second between 
V2,+ and vi + . We note that the bandgap size depends 
on the interparticle distance g inside the Wigner-Seitz 
cell but is independent on the number of cells M, and 
thus it is constant in the thermodynamic limit: in fact 
the quantization volume V oc l/y/M (in ID) gives that 
Q oc l/-\/M, so that the dependence on M in Eq. (|22|) 
cancels out. 

For g — > 0, in the limit of the monoperiodic array, 
one finds a single bandgap with size Auj = v\. + — v\ t ~. 
For g > (g < a) this interval is reduced: a frequency 
window opens inside the gap, where light is transmit- 
ted, and whose width is given by Auj g = ^ 2 .+ — u 2,-- 
The size of the two bandgaps is minimum at g = a/4, 
and it vanishes at this point when G\,q = Q2,Q- In 
this specific case, hence, the lattice becomes completely 
transparent. This is simply understood, considering that 
the bandgap results from an interference effect due to 
multiple scattering by all atomic planes, and it hence 
depends on the phase relations between the fields scat- 
tered by each plane. For this specific configuration, where 
g = A/4, 3A/4, the phase accumulated due to scattering 
of the first atom of the cell cancels out with the phase due 
to scattering by the second atom. As a result, the total 
phase accumulated from scattering with the two atoms 
of the cell is zero, and the medium hence behaves as it 
were completely transparent. Note that a monoperiodic 
array is also found for g = a/2, when all atoms have the 
same scattering properties. In this case the periodicity 
is halved, and the first BZ doubles, [—Go, Go]. The pho- 
tonic spectra of this system have been discussed in [8[ . 

These analytical results, obtained in a specific param- 
eter regimes, are confirmed by the results of the numeri- 
cal spectra, which are evaluated from Eqs. (|19[) by sum- 
ming over 40 BZs. The photonic spectra are shown in 
Figs. ElEUa), where the polariton dispersion relation is 
reported around q ~ for q = 0, 0.2a, 0.4a for several 
values of u>i, setting lo\ — u)2- The size of the bandgaps 
Alj + = i/2. + — and Aw_ = V2,- — as a function 
of g are displayed in Figs. [2!Ob), showing that the size 
of the gap is controlled by g, and it vanishes at Figures [4] 
and [5] display the photonic spectra when the atoms com- 
posing the Wigner-Seitz cell are of different species, in 
the case in which both interact with the probe but the 
resonance frequency of the respective dipolar transition 
is different. Figure 0] displays the photonic spectrum for 
the specific case in which one atomic transition is quasi- 
resonant, while the second is far detuned. In this case 
three photonic bandgaps appear, which vary largely as a 
function of g, and in such a way that while one is mini- 
mum, the other two are maximum, and vice versa. Fig- 
ure [S] displays the case in which the two atoms composing 
the cell are far detuned from the probe, with detunings 
of opposite signs. The spectrum is also characterized by 
three bandgaps. 



In this treatment we neglected atomic absorption, how- 
ever the evaluated bandgaps are significantly larger than 
the linewidth 7, and hence they can be experimentally 
observed. The effect of absorption is considered in the fol- 
lowing section in the framework of a scmiclassical model, 
where also finite size effects arc accounted for. 



IV. PROBE TRANSMISSION SPECTRA 

We now study the response of the biperiodic array of 
atoms to an external probe, and evaluate probe reflection 
and transmission as a function of the probe frequency 
considering the finite size effects. In this section we as- 
sume a three-dimensional array of atoms, where N lattice 
planes are aligned along the x direction. 

The coupling of an external probe to the atomic 
medium is described by the Maxwell equations for the 
electric field, in the presence of a dielectric in the region 
of space < x < L, with L = Ala. The equations for a 
probe at frequency oj p and wave vector k p = iv p /c, prop- 
agating along x, read 

[dl + k 2 p }E = -4nk 2 p P (23) 

where P(x, i) is the medium polarization. Continuity 
of the field and its spatial derivative at the interfaces 
imposes 

E(x = 0+) =E(x = 0_), 

d z F,(x = 0_) - d z E(x = 0+) = 4irk 2 P(x = 0) , ^ V 

and similarly at x — L. We consider the classical limit of 
these equations, which is found from the quantum model 
within an input-output formalism when the probe is a 
coherent state In particular, we calculate the reflec- 
tion and transmission coefficients of an incident field by 
means of the transfer matrix method. 

We consider a probe field propagating at normal inci- 
dence at x = to a single polarizable plane composed 
by atoms at resonance frequency ujj and linewidth jj . In 
the linear response regime, P = n s aj8(x)~E in Eqs. (|2"5]) 
and l|24p. where n s is the surface particle density and 

3 e 3 ( 25 i h l + i \ 

is the classical polarizability per particle, with Sj = 
ujj — ujv the detuning with respect to the atomic tran- 
sition [j. Denoting by E the amplitude of the inci- 
dent field, we write the reflected and transmitted fields 
as E r = rj~E e~ tkpX and E t = tj'Eoe lk '' x , where rj and tj 
are the transmission coefficients at a plane with polariz- 
ability a.j , and read 0] 

with £j = 2irk p n s aj. The transfer matrix Mj relates 
forward- and backward-traveling waves E^ , E^ on the 
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FIG. 2: (Color online) (a) Polariton dispersion for a ID bichromatic lattice for different values of g obtained by summing over 
40 BZ, when the atoms are trapped in a hollow fiber, whose fundamental mode is gaussian with waist w = 5/im 21]. The 
curves are evaluated for u>\ — u>2 = ^>o — IO7, considering the D2-line of 85 Rb atoms, with A = 780 nm and 7 = 2-7T x 6 MHz. 
The solid black line is plotted for comparison, and corresponds to the case in which atoms and field are not coupled (Qj = 0). 
(b) Photonic bandgap (numbered from lower to higher frequency) as a function of g/a, compared with the analytical prediction 
obtained from Eq. (f22|) . 
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FIG. 3: Same as in Fig. [2] but for ui = u)% — ujq ~ 5307- 



right-hand side of a plane to those on the left-hand side 



E+, E r , according to the relation 



E7 



(27) 



thereby automatically accounting for all interference ef- 
fects accumulated along the way. The transfer matrix 
for a single period is given by the product of the transfer 
matrix across the boundary of an atomic plane and free 
propagation for a distance dj 



1 







(28) 



For any transfer matrix M. = Y[j , generated by the 
product of single-period transfer matrices, the reflection 
and the transmission coefficients associated with the ma- 
trix elements are 



M12 
M22 



t 



(29) 



as is easily verified for the special case of matrix in Eq. 
([28]) . The transmission and reflection coefficients for the 
biperiodic lattice of M = N/2 planes is found as a func- 
tion of the transfer matrix = M-dx • M.d 2 f° r a 
single dimer, with d\ — g, di = a — g. In particular, the 
transmission coefficient reads 
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FIG. 4: Same as in Fig. [2] but for u)\ = luq — IO7 and cu2 — ^o + 5307- 
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FIG. 5: Same as in Fig. [5] but for ivi = o>o — 5307 an d W2 = wo + 5307- Note that gap 1 and 3 overlap. 



sin 8 



sin cos(nO) + i sin(n0)[£i£2 sin k p g — (1 — ^1^2) sin k p a — (£1 + £2) cos k p a] 



(30) 



(M n dld2 h2 

with O = acos[£i£2 cos k p g + (1 — £1^2) cos k p a — (£1 + £ 2 ) sin k p a] . The elementary cell dephasing can be written as 
function of the single slice dephasings Q\ 2 = acos(cos k p d\ 2 — £1,2 smk p dx 2) as following 



= acos 



cos(0! + 9 2 ) + sin 61 sin 9 2 — y sin 9i sin 6 2 — £1^2 sin k p g 



(31) 



showing that the term sin 2 k p g governs the difference between a simple addition of the dephasings Oj, as already 
pointed out in the calculation of the spectrum. 

If ImO > ReO, which occurs around lu p ~ wi,W2, in the limit n — > 00 the transmission coefficient (|30|) can be 
approximated by 

lim t - 2e-" Ime sine 

n->oo " sin9+£i£ 2 sinfcpg— (1 - £i£ 2 ) sinfc p a— (£1 + £ 2 ) cosfc p a 

I 



For vanishing values of Im0, the effect of multiple scat- 



tering cannot be factorized in a simple attenuation co- 
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FIG. 6: The absolute square of the transmission coefficient 
(left panel) and of the reflection coefficient (right panel) as a 
function of (lo p — lo)/j is plotted for a lattice of 85 Rb atoms 
with areal density n s = 5.7 ■ 10 _2 ^im -2 , trapped by laser 
beams detuned to the blue of resonance by IO7 and for 10 6 
atomic planes. 



efficient, and gives rise to interference structures as a 
function of g. 

In order to compare the transmission and reflection 
spectra with the results obtained for the photonic band 
structure, we consider a cloud of 85 Rb atoms in a quasi- 
one-dimensional geometry, confined in a very long lattice, 
taking N = 10 6 atomic slices. We fix the areal density at 
n s = 5.7 • 1CP 2 /im 2 , and we consider that the atoms are 
in the same internal state, so that L02 — and 72 = 71- 
In this section we will refer to optical lattices generated 
by laser beams with a frequency ujq = lui + IO7, where 
uj\ and 71 refer to the D2 atomic resonance at A = 780 
nm. Figure [5] displays the absolute square of the trans- 
mission coefficient (left panel) and of the reflection coef- 
ficient (right panel) of a probe beam traveling through a 
monoperiodic sequence of atomic planes (g = 0). For our 
choice of the parameters we find the edges of the gap at 
±4207. The width of the gap depends on our choice for 
n s . The spectrum of transmission agrees qualitatively, as 
a function of g, with the one shown in Fig. [3J where a 
strictly one-dimensional system was considered. Notice 
that the small peak in the transmission and the small dip 
in the reflectivity at the atomic frequency are a mark of 
scattering losses. 

In a bichromatic lattice, such as that shown in Fig. [T] 
the structure of the transmission and reflection spectrum 
is deeply modified by the multiperiodicity. Figures [7] and 
[5] refer to the cases g — 0.2a, and 0.24a respectively. 

In agreement with the results for the band structure at 
g — 0.2a in Fig. [2J at the center of the gap it appears a 
high transmission region divided in two symmetric parts 
by the absorption line (Fig. [7]). By increasing g (Fig. |5J), 
the width of the mini-band increases and the gaps become 
thinner. At g = 0.25a, the gaps close and the atomic lat- 
tice becomes transparent in this region of frequencies, ev- 
erywhere except for cu p ~ ljq, because of the absorption. 



FIG. 7: Same as Fig. H but g = 0.2 A. 



For this particular value of g, Md 2 — M-di — M. d , if 
£2 = £1 and neglecting absorption (the imaginary part of 
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FIG. 8: Same as Fig. H but for g = 0.24 A. 



V. BICHROMATIC LATTICE INSIDE A 
SINGLE-MODE CAVITY 

The observation of sufficiently large photonic bandgaps 
in free space requires a large number of lattice planes in 
a well controlled periodic structure, which is experimen- 
tally challenging. Nevertheless, observable effects can be 
found in small systems when coupling the atomic tran- 
sition, for instance, to the modes of a fiber [2l[ or of an 
optical resonator [231 ] . 

Let us assume that the dipolar transition of the atoms 
couples strongly with the single mode of a standing- 
wave optical resonator, which probes the system. The 
Hamiltonian for the dynamics inside the cavity is now 
H' = H dip + H c + H' int , where H dip is given in Eq. (J2J, 



Hr 



(33) 

describes the cavity field at frequency u> c , with a, an- 
nihilation and creation operators of a cavity photon, and 

~ gj cos(kxj + </?)er|a + H.c (34) 



3 



is the Jaynes-Cummings Hamiltonian, with k the cav- 
ity mode wave vector and gj the coupling strength, 
which depends on the dipolar moment Dj of the atom 
at position Xj. As a function of the cavity parameters, 
gj = y/q/(4wA)y/'y5tjj, with ? the scattering cross section 



in free space, A 



;/4 with w c the cavity mode waist, 



and Slo = 2itc/ L the free spectral range, with L the cav- 
ity length [l|| . The phase ip accounts for the dephasing 
between the cavity mode lattice and the atomic lattice. 
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The system is probed by an external weak coherent 
pump at intensity r\ and frequency ui p , which is coupled 
to the resonator. In this limit, we make the Holstein- 
Primakoff transformation and keep only the linear term. 
The resulting Heisenberg-Langevin equations of motion 
for cavity mode and spin wave operators read [24[ 



whether J\f is even or odd, respectively. The eigenfre- 
quencies of the homogeneous equations (for the specific 
case Si — 82 — A, 71 = 72 = 7) have the same form as 
in Eqs. (|36ti37[) . whereby now in Eq. (|3"Tj) the coefficient 
H reads 



11 = g\ cos 2 ip - 



g\ cos 2 (kg + cp). 



(39) 



-iS r a 



E 

Q£BZ,Q=k-G 

^ e -i(ke+ V ) dQ 



9i {e- iv b Q + 



<M 



.92 



b±Q = - 



d±Q 



(i<5i 
(162 
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2 



-VMgie 



±'up 



21 

2 



d±Q 



Mg 2 e 



We note that the eigenfrequencies in this case explicitly 
depend on g. This result is also found when the bichro- 
matic lattice is replaced by one single cell, by rescal- 
ing the coupling strength of each atom inside the cell 
, N as 9ea,j = vMgj. 

{60a,) now (jjscu^gg ^ ne intensity of the field at the cavity 

output as a function of the probe frequency u p for various 

JiBi, Q : parameters, by solving Eqs. (|35a[) - (|35c[) numerically. The 
^gg^quantity we study is the number of photons per unit time 
I{u! p ) — (o' ut a out ), where a out is the field at the cavity 
72^2, q j output, a ou t = \2~Ka — £, and the average is taken over 
the state of the system and the vacuum state of the e.m.- 
field outside the resonator [2J|. Hence, we find 



(35c) 



where S c = to c — lu p and k is the cavity lincwidth. The 
noise operators CW> &j,Qi have zero mean value and 
satisfy the relation <C(*)C(*') f > = (Bj.gWBj.Q^)') = 
5(t — t') (we assume the electromagnetic field in the vac- 
uum). 

Equations (|35a[) - (|35c[) describe the coupling between 
the cavity mode, at momentum fc, and the spin waves at 
quasi- momentum Q (inside the first Brillouin zone) , such 
that Q + G = k where G is a vector of the reciprocal 
lattice. We identify two relevant cases, when (i) k 7^ 
Nit /a, and (ii) k — Nit /a, where Af is an integer. 

For k j^z J\fir/a the system, composed by cavity poten- 
tial and Dichromatic lattice, is not periodic. The eigen- 
frequencies of the homogeneous equations can be simply 
found for the case <5i = 62 = A, 71 =72 =7, and read 



I(ojp) = 2n{a) a) 



2k V 2 (A 2 + 7 2 /4) 



(«7/2 - S C A + Mil) 2 + (Ak + 8 cl /2) 2 

(40) 

where we have used the steady state solution of Eq. (|35ap . 
In the strong coupling regime, when the cooperativity 
C ~ M7\L/2k7 > 1 [19j, the intensity I[u3 p ) at the cavity 
output exhibits two well defined maxima at the frequen- 
cies 



± 



Mil 



(41) 



v± 



A - i 7 /2 , 

5 C + A - i(« + 7/2) 



(36) 
(37) 



± 



A - ik + i- 
2 



Mil , 



where 



ll=(g 2 1 +g 2 2 )/2. 



(38) 



Here, the real part gives the position of the reso- 
nances, while the imaginary part gives the correspond- 
ing linewidth. The eigenmodes at frequency Vq are pure 
spin waves, and hence correspond to collective dipolar 
excitations which are decoupled from the cavity field. 
The eigenmodes at frequency v± are polariton excita- 
tions. We remark that the frequencies v± do not depend 
on g. 

For k = Af—, the system, composed by cavity poten- 
tial and bichromatic lattice, is periodic. The cavity mode 
couples to the spin waves Q — or Tv/a, depending on 



and cor resp onding to the vacuum Rabi splitting for this 
system 1!) . This can also be seen in Figs. H] and [TO] It is 
interesting to note that for A = and large cooperativity 
the cavity output field goes to zero as 1/C 2 . This is an 
interference effect, where the atomic polarization inside 
the cavity form a field equal and opposite to the driv- 
ing pump, such that the cavity field is effectively empty. 
Energy is in this case dissipated by the atoms. This be- 
haviour has been first predicted in 25j under the name 
" cavity induced transparency" . 

We evaluate the cavity transmission spectrum using 
the parameters of the setup in (26|, and consider 85 Rb 
atoms inside a resonator with length L = 85mm, fi- 
nesse T — 170000, loss rate k = 2ir x 21 KHz, beam 
waist w = 130 /im, and an average occupation per 
site n = 3000, that corresponds to an areal density 
n s ~ 5.7 • 10 _2 /zm -2 . The reflectivity of the cavity mir- 
rors is \r\ 2 ~ 1 — 1.8 x 10 -5 . The transmission spectrum 
is calculated assuming that N — 200 planes (M = 100) 
are confined inside the resonator. The positions of the 
peaks correspond to those predicted by Eq. (|4Tj) by tak- 
ing into account multiple occupancy of the lattice sites 
rescaling the coupling strengths as gj — > v^5i ■ Figure [5] 
displays the squared transmission as a function of the 
probe frequency for the case if = tt/2 and the values 
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FIG. 9: (Color online) The absolute square of the trans- 
mission coefficient as a function of (w p — u>o) (in units of 
7) is plotted for a lattice of 85 Rb atoms with areal density 
n s = 5.7 ■ 1CF 2 /im~ 2 , trapped by laser beams detuned to the 
blue of resonance by 10 7, for 200 planes, in the presence of 
a cavity of length L — 85 mm and finesse T = 170000, and 
for <p — n/2. Note that for g = the atoms are trapped at 
the nodes of the resonator and hence do not interact with the 
cavity field. 
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FIG. 10: Same as Fig. but for tp = 0. 



inside the primitive Wigner-Seitz cell. We have derived 
a model describing light propagation for a weak probe, 
and its response to probe propagation in free space and 
inside of a cavity. We have found that, depending on g, 
in free space the system may or may not exhibit pho- 
tonic bandgaps about the atomic frequency. This is a 
peculiar property, which makes the biperiodic crystal dif- 
ferent from the monoperiodic one, always exhibiting a 
bandgap at the atomic frequency. In case there are pho- 
tonic bandgaps around this value, they occur in two or 
more ranges of frequencies. For a finite crystal, relevant 
effects can be observed when the atoms are confined in- 
side an optical resonator. Here, the interparticle distance 
g inside the primitive Wigner-Seitz cell determines the 
properties of the transmission spectrum of a probe at 
the cavity output. 

Our study is based on a full quantum model for the 
light. In this paper we have focused on the intensity 
of the transmitted and reflected light, in the future we 
will study higher order coherence of the scattered light. 
On the basis of studies made with two atoms inside 
a cavity [27| . we expect that, when considering satu- 
ration effects, the biperiodic optical lattice can act as 
nonlinear-optical medium, whose properties may be con- 
trolled by the interparticle distance g. In the linear re- 
sponse regime, it can be interesting to study higher order 
coherence of the scattered light for various states of mat- 
ter inside the potential, with the aim of determining the 
quantum state of the matter [l(|- This may allow, in 
particular, to detect experimentally novel states of mat- 
ter realized in bichromatic optical lattices [H, [TH, HU • 



g = 0,0.2a, 0.4a. Note that the case g = corresponds 
to all atoms at the antinodes of the resonator, and it is 
hence equivalent to the situation in which the cavity is 
empty. Figure [10] displays the transmission spectrum is 
the optical lattice trapping the atom is shifted so that 
ip = 0, showing that the form changes substantially. By 
varying ip, hence, information on the interparticle dis- 
tance in the Wigner-Seitz cell can be gained. The min- 
imum at ujp = u>o, corresponding to the cavity induced 
transparency behaviour, is here visible. 

VI. CONCLUSION 

The photonic properties of biperiodic optical lattices 
are critically determined by the interparticle distance g 
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